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The magnetic field distribution in the barrier of small planar Josephson tunnel 



Q I junctions is numerically simulated in the case when an external magnetic field is 



applied perpendicular to the barrier plane. The simulations allow for heuristic an- 
alytical solutions for the Josephson static phase profile from which the dependence 
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f^ I of the maximum Josephson current on the applied field amplitude is derived. The 



most common geometrical configurations are considered and, when possible, the 
theoretical findings are compared with the experimental data. 
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I. INTRODUCTION 



The static (and dynamic) properties of a planar Josephson Tunnel Junction ( JTJ) are well 



y. On 



understood when an external magnetic field is uniformly applied in the junction plane 

the contrary, very little is known when a uniform magnetic field is applied perpendicularly 

to the barrier plane. The main reason why, since the discovery of the Josephson effect in 



1962, only few papers have dealt with a transverse magnetic field|2| 



3||^, is due to the fact 



that demagnetization effects imposed by the electrodes geometry are awkward to take into 
account. In a recent paper [5| we provided an experimental proof that a transverse magnetic 
field can be much more capable than an in-plane one to modulate the critical current Ic of a 
planar JTJ with proper barrier and electrodes geometry requirements. It is also possible to 
design the JTJ electrode geometry in such a way that it is totally insensitive to a transverse 
field. The possibility to have on the same chip JTJs having different sensitivities to an 
externally applied field can be very attractive in practical applications. 
In this paper we push our analysis further by resorting to numerical magnetostatic simu- 
lations to find the field distribution H in the barrier plane of those JTJs having the most 
common rectangular and annular geometries. Once H is found empirically, the Josephson 
phase d), which is the difference between the complex wavefunction phases in the electrodes, 

n 

can be obtained from the Josephson equation[6|]: 

V0 = — H X n, (1) 

<i>o 

where n is a unit vector normal to the insulating barrier separating the two superconducting 
electrodes, /xq is the vacuum permeability and $o = h/2e is the magnetic flux quantum. If 
the two superconducting fllms have thicknesses (ii 2 and London penetration depths \li,2 
and tj is the barrier thickness, then the effective magnetic penetration d^. is given 



di ^ d: 



2 



de = tj + Ali tanh — h \l2 tanh , 

•^Ali /Al2 

which, in the case of thick superconducting fllms {di » Xu), reduces to d^ ~ Ali + Al2 

(being always di » tj). 

In Cartesian coordinates, assuming that the tunnel barrier lies in the x — y plan, Eq{T] reduces 

to: 



d(p{x,y) d(l){x,y) , . 

^^«-i7„ -^^ocH^. (2) 

For a planar JTJ with a uniform Josephson current density Jc whose dimensions are smaller 
than the Josephson penetration depth Xj = ^Jh/2e^QdeJc-i the self-induced field associated 
with the bias current can be neglected and the Josephson phase must satisfy the two- 
dimensional Laplacian equation |6|: 

^ + ^ = 0, (3) 

dx^ dy"^ 

with proper boundary conditions related to the value of the magnetic field components i^a; 
and Hy on the junction perimeter. It was first pointed out in 1975 3| that in a transverse 
applied field H = H_lz, the in-plane components H^ and Hy are ascribed to surface de- 
magnetizing currents js = z x H feeding the interior of the junction. Since these currents 
mainly flow on the film edges, the largest sensitivity to a transverse field occurs when the 
junction is formed at the film edges. On the contrary, if the barrier is placed well inside 
the superconducting films, the effect of a transverse field vanishes. Our task consists of 
numerically evaluating the field line distribution in the barrier plane, from which we deter- 
mine an empirical analytical expression (j){x, y) for the phase profile which satisfies EqOl 
Such a phase profile will allow the computation of the transverse magnetic diffraction pat- 
tern Ic{Hj_) for small JTJs having different geometries and to compare it with experimental 
data, if available. This is achieved by recalling that the maximum Josephson current is: 

Ic = Jo a/ (sin 0) 2 -f (cos 0)2, 

in which the brackets () denote spatial averages over the junction area. Throughout the 
paper we assume that the applied transverse field is everywhere much smaller than the 
critical field which would force the films into the intermediate or normal state, i.e., that the 
superconductors are always in the fiux-free Meissner regime. 

II. MAGNETOSTATIC SIMULATIONS 

In general, magnetostatic problems are based on the magnetic vector potential. However, 
where no electrical currents are present, the problem can be conveniently solved using the 



scalar magnetic potential. In fact, in a current free region V x H = allows the introduction 
of a scalar potential Vm such that H = — VVm- Using the constitutive relation B = yUoA^rH, 
we can rewrite Maxwell's equation V • B = in terms of Vm'- 

- V ■ (/io/^r VV„) = 0, (4) 

in which the magnetic relative permittivity fir is spatially dependent. We assumed that the 
superconducting electrodes are thicker than their London penetration depths (rfj >> Xu), 
so that the London equation reduces to B = everywhere inside the superconductors, i.e., 
Hr = (perfect diamagnetism) and the normal component of the magnetic flux density 
vanishes at the boundary (n ■ B = 0). In the opposite limit, the films would become 
transparent to the transverse field and, in turn, the junction would lose its sensitivity to 
the transverse field. A uniform applied magnetic field H = H^z is taken into account by 
imposing that sufficiently far away from the junction is Vm = —Hj_z. All the simulations 
presented in this paper were carried out setting H± = lA/m. 

As a consequence of the definitions of Eql2], it is straightforward to show that Eql3] requires 
that d'^Vm/dxdy = d'^Vm/dydx. Further, more importantly, we have: 

0(x, y)(x j dx-^ = - J dy-^. (5) 

The numerical solution of Eqj4] was implemented in the COMSOL Multiphysics 3D Elec- 
tromagnetics module for JTJs having different rectangular and annular geometries. Models 
with large geometric scale variations are always problematic to mesh, in particular if they 
contain thin layers with large aspect ratio. Therefore, one caveat of our modeling is that, in 
order to keep the number of mesh elements within the PC memory handling capability, the 
separation between the superconducting films, i.e. the tunnel barrier thickness tj, could not 
be set to realistic values for a Josephson tunnel barrier 0{lnm). Our numerical modeling 
was tested against the magnetic field distribution around a superconducting disk (with ra- 
dius R larger than its thickness d) in the plane z = 0, centered on the z axis and immersed in 
a field H = Hj^z. More precisely, the radial dependence of the tangential field Ht on the disk 
surface followed to a high accuracy the well known expression |8|] Ht{r) oc H^r/\/E? — r^ 
everywhere except at the disk border, where the inverse square root singularity was replaced 



by a finite value H proportional to the square root of the disk aspect ratio H = H±^ ^jRjd 



Q. 




FIG. 1: Sketch of a overlap geometry junction . The center of the junction coincides with origin 
of our coordinate system. 

This example is indicative of the fact that, in general, the magnetostatic response of any 
superconducting film structure is markedly dependent on the film aspect ratio. 

III. RECTANGULAR JUNCTIONS 

A. Overlap type junctions 

We begin our analysis with a JT J obtained by the superposition of the extremities of two long 
and narrow parallel superconducting electrodes with equal widths. This so-called overlap 
geometry is depicted in FigH] for a square junction, i.e. W = L. The tunnel barrier lies 
in the z = plane and its center coincides with the axis origin. Further, it has a length 
2L along the x-direction and a width 2W along the y-direction. In the simulations the 
electrodes have a thickness d = lOfim and are tj = Ifim apart. The film width 2W and 
the overlapping length 2L where varied in order to treat barriers with different aspect ratios 
/3 = L/W. Figsl2^-c show the numerically obtained Vm solutions in the barrier area of three 
overlap junctions having the same width 2W = SOfim, but different lengths 2L = 20, 80, and 
320fim. By analyzing the properties of such plots we aim to infer an empirical, physically 
acceptable analytical form for Vm{x,y,z = 0). We observe that, for any value of /?, the 
scalar potential in the barrier is symmetric with respect to the x-axis and antisymmetric 
with respect to the y-axis. In other words, the expression Vm{x, y) we are looking for has to 
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FIG. 2: (Color online) Numerically obtained magnetic scalar potential Vm (in A) for three overlap 
planar Josephson tunnel junctions having the same width 2W = 80fj,m, but different lengths: a) 
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FIG. 3: (Color online) Numerically obtained in-plane magnetic field components in A/m for a 
square overlap junction having 2W = 2L = 100/um in a transverse externally applied field H± = 
lA/m. Color plots for a) Hx{x,y,z = 0) and b) Hy{x,y,z = 0). 

be an odd function of x and a even function of y. Furthermore, we note that the potential 
decays from the junction corners over a distance W, being mostly null when L > W (or 
/3 > 1). We have checked that the following ansatz: 



Vm{x,y) = WH cos qy 



sinh qPx 



(6) 



sinh/3 

in which g is a fitting parameter near unity, allowed us to reproduces the plots of Figs|2^-c 
at a better than qualitative level. In fact, for q = 1 the relative difference between the 
simulation output and the proposed expression being everywhere less than ±15% and the 
q value that minimized the error was q ~ 0.9. We have introduced the normalized units 
X = x/L and y = y/W (note that /3x = x/W). In the last equation, again H oc H^^/wJd, 
with a proportionality constant of order of unity which slightly increases when the barrier 
thickness tj decreases. Unfortunately, recalling the comments of the previous section, we 
cannot be more precise on this point. 

Now we focus our attention on the components of the magnetic field in the barrier plane 
Hx{x,y) = —dVm/dx\z=o and Hy{x,y) = —dVm/dy\z=o {H^ being identically null all over 
the barrier area). They are shown in FigJSti and b, respectively, for the particular case 
W = L = lOOfim. We like to specify, at this point, that the same plots obtained from 
numerical simulations based on the vector, rather than scalar, potential differed by no more 



8 

than ±10%, the discrepancy being larger at the barrier edges. From Eq|6]with q = 1, the 
following analytical expressions are derived: 

TT /- -\ TT ^ cosh l3x ,_, 

H^ix,y) = -Hcosy—^, (7) 

H,ix,y) = Hsmy—^. (8) 

The physical meaning of the last expressions is that for (3 = L/W > 1, the magnetic field 
lines are confined to the corners of the junctions at a distance W and most of the field lines 
entering the junction at x = ±L are bent by 90° and leave at y = ±W. In the opposite 
limit, cosh/?x ~ 1, so the x-dependence of H^ disappears, meaning that all the field lines 
entering the barrier at, say, x = —L exit at x = L (or viceversa). Further, we notice that 
while the x-component is negative all over the barrier area, the y-component symmetrically 
spans from negative to positive values. Due to the linearity of EqH] e the system symmetry 
with respect to the z = plane, if the direction of the transverse field is reverted, then 
Hx and Hy simply invert their sign. The magnetic field line distributions in the junction 
barrier corresponding to the scalar potentials of FigsE^i-c are shown in FigsH^-c. Similar 
plots based on the previous analytical expressions would be practically undistinguishable 
at the picture resolution level, therefore, they will not be shown. From the magnetic field 
distributions we expect that, for a given junction area LxW, the critical current Ic of a 
planar JTJ with pure overlap geometry [L < W) modulates much faster than that of a 
sample with pure in-line geometry {L > W). At a first sight, it might seem that the effect 
of a transverse field is qualitatively similar to that of an in-plane field applied along the film 
direction, i.e. along the x-axis, in our case. However, this is not true at a quantitative level 
because, in general, Hx{x, ±W) is not constant in a transverse field. 

Inserting Eq|6] in any of the Eqj5], we derive an approximate analytical expression for the 
Josephson phase profile: 

,/- -\ „cosh/3x , . 

0(x,y) = /ismy— ^, (9) 

where h = 2'n'deWfioH/^Q is a dimensionless parameter proportional to the applied trans- 
verse field amplitude H± through H. It is easy to verify that the last expression, in which 
we have omitted an integration constant (pQ, satisfies Eq|3l 
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FIG. 4: Magnetic field lines inside three overlap planar Josephson tunnel junctions having the same 
width 2W = 80/im, but different lengths: a) 2L = 20/im (/? = 0.25), b) 2L = SO^um (/? = 1) and 
c) 2L = 320//m (/3 = 4). 

With an odd function of y, then (sin0) = 0; therefore, the magnetic pattern /c(/i) reduces 
to: 



/,(/i) = /,(0) I dx I dycos(l){x,y). (10) 

JO io 

FigjS^-c show the computed transverse magnetic patterns for the three values of the barrier 
aspect ratio /? used before (4, 1 and 0.25). As expected the response to a transverse field 
is very weak for an in-line JTJ, the first minimum occurring at /i ~ 74 for (3 = 4. The 



10 



secondary pattern maxima become more pronounced for a pure overlap geometry. However, 
in the limit L « W, all the above equations lose their validity when the overlapping length 
becomes comparable with the film thickness. 
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FIG. 5: Computed transverse magnetic patterns Ic{h) for an overlap junction with different L/W 
ratios: a) inline junction L = AW, b) square overlap L = W and c) pure overlap junction W = AL. 



It is important to stress here that we are dealing with electrically small JTJs, therefore 
the diff'erent shapes of the transverse magnetic pattern is a direct consequence of the dif- 
ferent distribution of the surface screening currents (and not of the applied bias current). 
Unfortunately, there are no data available in the literature to check the validity of our the- 
oretical magnetic diffraction patterns for a small overlap JTJ formed by films having the 
same widths. In fact, the experiments reported by Rosestein and Chen in 1975 ^ refer to 
an overlap JTJ formed by two thick Ph electrodes of unequal widths {2W = 0.74mm and 
2W' = 1.00mm) and a common overlay region of 2L = 0.35mm. It is quite evident that, for 
such geometrical film configuration, the symmetry with respect to the y-axis is broken and 
Eqj9] is unable to correctly describe the magnetic field (and screening currents) distribution. 
FigEK and b show, respectively, the result of numerical simulations carried out for the spe- 
cific electrode configuration of Ref. [2| and the corresponding Ic{h). According to Ref. [3|, we 
believe that difference between the experimental data of Ref. 2[ and the numerical prediction 
of Ref. p] valid only for the specific case P = 0.5, arises from the unequal widths of the films 
in the experiment. Indeed, the magnetic diffraction pattern reported in Ref.[3j is of a piece 
with the curve in Figl6b. 

We conclude this section considering that, for unidimensional overlap junctions for which 
W > Xj > L, being ^ << g-|, then the Josephson phase has to obey to the equation first 
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(a) (b) 

FIG. 6: (Color online) a) Numerically obtained scalar potential for an overlap junction obtained 
by the superposition of two films of unequal widths 2W = 0.74:mm and 2W' = l.OOmm. The 
overlapping distance is 2L = 0.35mm, as for the sample quoted in Ref.[2j. b)Computed magnetic 
pattern Ic{h) for an overlap junction with aspect ratio L/W = 0.5. 

introduced by Owen and Sacalapino[lO|: 



rf^0 1 . , 



(11) 



when an in-plane external field is applied along the x-direction. In fact, FigJT^ shows the 
comparison between the diffraction patterns measured in a parallel and transverse field of a 
Nb/Alox/Nb overlap-type junction with Xj ~ 50fim whose length is 500fim, while the width 
is equal to Afim. The base and top electrode widths are 540 and 506^m, respectively. The 
junction geometry is depicted in FigJ7|D. We observe that the two experimental datasets 
almost overlap, when a factor scale of about 40 is applied on the abscissae, meaning that 
the sample is much more sensitive to a transverse field rather than an in-plane one. 



B. Cross type junctions 

Cross geometry JT Js are formed by the superposition of two perpendicular superconducting 
electrodes, as depicted in FiglHl The static properties of such junctions in a transverse 
magnetic field were analyzed by Miller et al. J], but only in the particular case of equal 
film widths 2L = 2W . They proposed, as an approximate solution of EqjSl a phase profile 
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FIG. 7: a) (Color online) Comparison between the diffraction patterns measured in a parallel 
(black bottom axis) and transverse (blue top axis) field of a Nb/Alox/Nb overlap- type junction 
with Xj ~ SO/um whose length is SOO/um, while the width is equal to 4:fj,m. b)Geometry details: 
the base (red) and top (blue) electrode widths are 540 and 506fj,m, respectively. The barrier area 
is delimited by the black rectangle. 




FIG. 8: Sketch of a square cross type junction. The center of the junction coincides with the axis 
origin. 

0(x, y) oc xy (corresponding to a saddle shaped scalar magnetic potential Vm{x,y) oc x^ — 
y'^ and to a monotonically decreasing Ic{H±).) We want to generalized these results for 
junctions with non unitary aspect ratios /3 = L/W. FigsIHti'-c display the numerical solutions 
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FIG. 9: (Color online) Numerically obtained magnetic scalar potential Vm (in A) in the barrier 
plane for three cross type junctions having the same width 2W = lOO/xm, but different lengths: a) 
2L = 100/im {(3 = 1), b) 2L = 40/xm (/3 = 0.4), and c) 2L = 20fim (/3 = 0.2). The external applied 
field is HjL = lA/m. 



of EqJUfor three cross junctions having the same width 2W = 100 fim, but different lengths 
2L = 100, 40, and 20fim. We observe that, for any value of f3, the scalar potential in 
the barrier is four- fold symmetric meaning that the empirical expression Vm{x,y) we are 
looking for has to be an a even function of both x and y. Further, Vm always vanishes at 
the junction corners and, as the junction length shrinks, the scalar potential distribution 
gets more and more uniform over the barrier area. A careful analysis of the scalar potential 
plots in Figs|9K-c, led us to the following expression: 



T/ f- -\ ir /TUT f .coshg/3x ^ cosh qy/(3 \ 



(12) 



14 

in terms of normalized variables. Here again g is a fitting that can be comfortably set equal 
to unity. More specifically, with q = I, the relative difference between the simulation output 
and the heuristic expression of Eqll2] was numerically found to be everywhere less than 
±15% although it was minimized by g ~ 1.2. The proposed expression is made up by two 
terms which can be seen as the contributions from the two electrodes. When L = W, the 
two terms have the same weights (sechl ~ 0.65), but, for, say, P < I, the weight of the first 
term is larger than the one of the second term and viceversa. Further, in the limit (3 << 1 
the first weight saturates to unity while the second vanishes. 
From EqlT2]with unitary q, the Josephson phase profile can be easily derived: 

,.. .^ , /. .sinh/3x . „sinhy//?\ 
0(x, ,) = h [sm ,— ^ + sm x—^^) , (13) 

where h = 2TTdeVWLfioH/^Q and with Eq|3] being identically satisfied. We begin with the 
observation that setting (3 = 1 and retaining the first two terms in the Taylor expansion of the 
trigonometric and hyperbolic functions, EqsilS] and [12] reduces to 0(x, y) oc xy — x^y^/36 
and Vm{x,y) oc x^ — y^, as it should be. Further, upon the inversion of /?, (p{x,y) = 
(j){y,x), meaning that the solutions for two junctions having reciprocal aspect ratios differ 
by a rotation of ±90°. FigsITOk-c show the magnetic field distributions in the barrier area 
corresponding to the scalar potentials shown in Figsj9K-c. 

Again, with (sin0) = 0, the magnetic diffraction pattern for a cross junction in a transverse 
magnetic field are found on inserting the expression above in EqJTOl FigJTT] shows Ic{h) for 
the three values of the barrier aspect ratios considered in Figsj9ti-c, i.e., (3=1, 0.4, and 0.2. 
For the considerations above, the red and black curves in Fig|TT]also represent the Ic{h) for 
(3 = 2.5 and 5, respectively. We come to the interesting result that for cross junctions in 
a transverse field the critical current decreases monotonically with the field amplitude H± 
and, for large fields {h » 1), Ic{h) oc 1/H± (see the log-log plot in the inset of Fig JTTj) . The 
experimental transverse pattern presented in Ref.jJ] bears strong resemblance to the Idh) 
obtained for (3 = 1. 

IV. ANNULAR JUNCTIONS 

In this section we will examine the behavior of small annular JTJs in the presence of a 
transverse magnetic field. Denoting the inner and outer ring radii, respectively, as Tj and 
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(b) (c) 

FIG. 10: Numerically obtained magnetic field distribution in the barrier area of three cross-type 
planar Josephson tunnel junctions having the same width 2W = WOfim, but different lengths: a) 
2L = lOOfim (/3 = 1), b) 2L = 40^m (/? = 0.4) and c) 2L = 20fim (/? = 0.2). 



To, we assume that the annular junction is unidimensional, i.e., the ring mean radius r 
{ti + ro)/2 is much larger than the ring width Ar = {vo — rj)/2. 
Using polar coordinates, the Josephson equation EqUJcan be split into: 
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FIG. 11: (Color online) Computed transverse magnetic patterns Ic{h) for a cross junction with 
different L/W ratios /3 = 1, /3 = 0.4 and /? = 0.2. In the inset the log-log plot shows that for large 
fields Icih) oc 1/h. 




(a) (b) 

FIG. 12: (Color online) Sketches of the two types of annular JTJs considered in this paper, a) 
Lynghy type geometry made by two films whose widths match the ring diameter, b) Asymmetric 
annular junctions made by two films with unequal widths. 



kHb 



dcj) 



KjII f . 



(14) 



dr rdO 

where Hr and He are the radial and tangential components of the magnetic field in the ring 
plane, respectively and k depends on the electrodes geometrical configuration |l4|. With 
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the annulus unidimensional, we can neglect the radial dependence of the Josephson phase 

0(r, 6) = (j){r, 6) and, henceforth: 



= -nr / dOHrir, 6) + const. (15) 

In the well known case of a spatially homogeneous in-plane field Hy applied in the direction 
of ^ = 0, then Hr = H\\ cos 6 (and Hr = H\\ sin^), so that the last integral yields 12| |: 



= /isine + 0o, (16) 

where h oc H\\ and 0o is an integration constant. Assuming that the Josephson current 
density Jc is constant over the ring circumference, the Josephson current through the barrier 
is obtained by: 

uh) = ^-^ r dosiTim- 

in which /(O) = Jc27rfAr is the maximum junction critical current which occurs in zero 
field. As far as 0(6') is an odd function (when 0o = 0), the calculation of the maximum 
critical current reduces to the following integration: 

^h) = ^^ r decos<j){e). (17) 



Inserting given by EqJT6| we obtain for the maximum critical current 14] 



/,(/l)=/e(0)|jo(/l)|, (18) 

in which Jq is the zero order Bessel function (of first kind). The periodic conditions for the 
phase difference (p and its angular derivative around an annular junction are: 

+ 27r) = (p{9) + 2™, (19) 

+ 2vr) _ dm 

de ~ do ' ^^^^ 

where n is an integer corresponding to the net number of fluxons (i.e., number of flux- 
ons minus number of antifluxons) trapped in the junction at the time of the normal-to- 
superconducting transition. Eqs. flT^ and fl^Ul) state that observable quantities such as the 
Josephson current (through sin0) and the radial magnetic field (through d(f>/d9) must be 
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single valued upon a round trip; they were derived in Ref. IJ] starting from the fluxoid 

quantization. 

EqUniand EqIT8]hold under the assumption that there are no fluxons trapped in the barrier; 

however, they can be easily generalized to the case of n 7^ trapped fluxons. In such case, 

EqJT6] changes to: 



= hsme + ne + (j)o, (21) 

in which the linear term in Eql2T] takes into account the phase twist due to the presence of 
the trapped fluxons, being that the ring circumference is smaller or comparable to the fluxon 
rest length. Carrying out the integration in EqUT] with given by Eqj2T] and maximizing 
with respect to 0o, we get: 

i:{h) = lM\Uh)\, (22) 

in which J^ is the n-th order Bessel function. EqsiTS] and [22] have been experimentally 
verified in a number of papers. 

n 

A Lyngby type annular JTJ, firstly reported in 1985 by Davidson et a/.[15i], is obtained by 
two films having the same width, as schematically depicted in FigJT2k. Further, FigJT2b 
shows a different kind of annular JTJ for which the film widths are quite different: we will 
call it asymmetric annular junction. At the end of this section we will present experimental 
data for such asymmetric geometrical configuration. We have carried out magnetostatic 
simulation for the two annular geometries depicted in FigsHSl when the applied field is 
transverse. Only the case of no trapped fluxons was considered, corresponding to zero net 
magnetic flux through the superconducting holes. Contrary to the case of the rectangular 
bidimensional JTJs considered previously, now we do not have to know the magnetic field 
distribution in the junction plane, but, by virtue of EqJT^. we can limit our interest to 
just the angular dependence of the radial magnetic field Hr{r, 6). In our simulations we set 
r, = 40/xm and Tq = 50fim, so that r = 45/im. For the asymmetric configuration the film 
widths were chosen to be 2W = lOOfim and 2W' = 200fim. 

Postprocessing the simulation outputs we found out that, in the case of Lyngby geometry, 
Hr follows very closely a sinusoidal dependence on 6, as shown in FigsJTSh: more specifically, 
by choosing the angle origin in such a way that 6 = corresponds to the positive x-axis 
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FIG. 13: (Color online) Angular dependencies of the radial magnetic field Hr{6) for a small annular 
junction in a transverse field having: a) Lynghy type geometry sketched in Figll2b. and b) asym- 
metric geometry sketched in Fig|12b. In the numerical simulations the amplitude of the transverse 
field was set to lA/m. 

direction, we have H^ oc cos^, exactly as if the magnetic field were applied in the ring plane. 
By integration we get EqEU again with h depending on the geometrical film configuration 
and being proportional to the transverse field amplitude iJ^. We come to the remarkable 
conclusion that the diffraction pattern of an electrically small annular junction with no 
trapped flux in a transverse field follows the zero order Bessel function behavior, as if the 
field were applied in the barrier plane. 

The situation is quite different when we consider asymmetric annular junctions. In fact, as 
shown in FigJTSb. it is quite evident that now the slope of the radial field changes abruptly 
for 9 ~ 7r/2 and 9 ~ 37r/2 resulting in a periodic asymmetric ratchet-like potential dHr/d9. 
We have numerically checked that to a high accuracy j^ d9Hj.{9) = 0, as it should be when 
no fluxons are trapped in the junction. In order to correctly reproduce Hr{9), we have to 
consider higher 9 harmonics. It was found that a truncated Fourier expansion cast in the 
form: 



i7^(^) oc cos + 27 cos 20 + 35 cos 3^ (23) 

can satisfactorily fit our numerical findings. The two fitting parameters 7 and S can be 
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ascribed to two degrees of freedom in the layout geometry: the ratio of the top and bottom 
film widths and the distance from the junction to the edge of the bottom film. Eqj23] with 
7 = 0.11 and 6 = 0.085 is shown as a solid red line in FigJT3b. 

By integrating EqJT^ with Hr given by Eqj23l we get an approximate expression for the 
angular phase dependence: 



:i h{sme + -fsm2e + 6sm3e) (24) 

in which still (j){-0) = -<f){6), but the symmetries 0(7r/2-6') = (j){n/2 + e) and (j){3n/2-e) = 
0(37r/2 + 6) are now lost. Since (f){6) is an odd function, Eq{T7] allows us to calculate the 
magnetic diffraction patterns corresponding to the above non-sinusoidal phase profile, even 
in the case when a term n6 is added to account for the presence of n trapped fluxons. 
It turned out that, while for n = the effects of the 7 and S terms tend to cancel each 
other, resulting in a zero-order Bessel function behavior as in EqJTSl for n 7^ we found a 
marked departure from the n-th order Bessel function dependence of Eqj22l These results 
are supported by experimental results for an asymmetric annular junction (r = 80 fim and 
Ar = 4yum) made by unequal width films: the base electrode width is 540/im and the top 
electrode width is 170 fim. For such layout, the numerical analysis of the angular radial field 
dependence yielded the best fit values 7 = 0.19 and 6 = 0.078. In the FigsHHi and b we 
show, respectively, the experimental diffraction patterns (dots) for such junction without 
trapped fluxons and with n = 1 trapped fiuxon. The experimental data can be fitted very 
nicely by the theoretical expectations (solid red lines) obtained inserting the above 7 and S 
values in Eqj24l 

We observe that when no fluxons are trapped in the asymmetric annular junctions the 
transverse pattern is definitely symmetric with respect to the inversion of field direction and 
is barely distinguishable from the pure Bessel one; further, we stress that the same sample 
measured with an in-plane field applied in the ^ = direction showed again a Bessel like 
pattern, but the response to the applied field was about 25 times weaker. 
On the contrary, with n = 1 the transverse magnetic diffraction pattern loses its symmetry 
with respect to the field amplitude, i.e. Il{—h) 7^ Ilih). Furthermore, both in the exper- 
iments and in the calculations, it turns out that I^{~h) = I~^{h); in other words, if we 
invert both the field and fiuxon polarities we obtain the same magnetic diffraction pattern. 
This result was obtained and exploited in the context of a detailed investigation of the sym- 
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FIG. 14: (Color online) Transverse diffraction patterns for an asymmetric annular junction. The 
dots are the experimental data, while the sohd red hues are the theoretical expectations obtained 
inserting 7 = 0.18 and 5 = 0.08 in EqlJH a) no trapped fluxon; b) one trapped fluxon. 

metry breaking during fast normal-to-superconducting phase transitions of annular JTJs 
recently published 111] . Among other things, it has been experimentally and theoretically 
demonstrated that when a small transverse field is applied to the ring during the thermal 
quench the probability to trap a Josephson fluxon can be very close to unity, the fluxon 
polarity depending on the field polarity. The ability to easily discriminate between a fluxon 
and an antifiuxon can be conveniently exploited in the recently proposed Josephson-vortex 
qubits experiments with ring and heart-shaped JTJs 16] . The asymmetry of the magnetic 
diffraction pattern can be very simply ascribed to the ratchet-like pot ential whose effect on 
the fluxon dynamic properties has been fully investigated recently |17l]. 
We conclude this section by remarking that the angular dependence of both the radial and 
tangential magnetic field components in the barrier of a annular JTJ do not change if the 
circular hole is removed from one of the electrodes. This is supported by both numerical 
simulations and experimental data[l8|. Indeed, when the ring shaped barrier is formed 
between a holed film and a singly connected one, the Josephson fluxon polarity is univocally 
related to the polarity of the quantized flux threading the hole. 



22 



V. CONCLUDING REMARKS 

The transverse magnetic patterns of electrically small Josephson tunnel junctions have been 
derived numerically by solving the magnetostatic problem for different geometrical configu- 
rations of the junction electrodes and of the barrier. More specifically, from the numerical 
analysis of the magnetic scalar potential produced in the barrier plane by the demagnetizing 
currents circulating on the electrode surfaces we derived approximate and simple expres- 
sions for the Josephson phase distribution in the barrier area, which, in turn, permitted to 
calculate the junction critical current. Such calculations show, among other things, that 
for rectangular barriers the modulation of the maximum critical current never follows the 
Fraunhofer behavior typical of a field applied in the barrier plane; further, Ic{H±) strongly 
depends on the barrier aspect ratio L/W. On the contrary, the critical current modula- 
tion in a transverse field of annular JTJs without trapped fluxons is fairly close to the one 
corresponding to a parallel field, although it can be much faster when the field is perpen- 
dicular. When the film configuration of the annular junction is asymmetric, then the static 
properties depend on the polarity of the transverse field and of the trapped fluxons. It's 
worthy to mention that our calculations were carried out assuming that the junctions were 
not biased. However, in order to measure the magnetic diffraction patterns one needs to 
supply a transport current by an external source. As far as the JTJ is electrically small, as 
in cases considered in this paper, the effect of a non uniform current distribution through 
the barrier is negligible 1|] . Nevertheless, to exclude flux from the electrodes interiors, a 
self- field that wraps around the films is generated [19!] whose effect on the Josephson phase 
distribution is largest when the current is largest. This situation typically occurs when the 
applied field is small (or absent) regardless of its orientation with respect to the barrier 
plane. As the external field amplitude grows, the relative effect of bias induced screening 
currents decreases, and disappears when the field amplitude is such that the critical current 
is zero. 

It is important to stress that the static properties of a small JTJ in a transverse field is 
strictly related to the film layout. In the case of junctions formed in a windows between two 
films which completely overlap each other near the junction itself, the circulating currents on 
the film interior surfaces are symmetric with respect to the barrier plane and result in a zero 
magnetic field; consequently such JTJs will remain totally insensitive to a transverse field: 
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this holds for overlap type and annular geometry JTJs shown, respectively, in Fig{T] and 
FigJT2k when one of the electrodes is rotated by 180°. As mentioned in the Introduction, 
we also expect a very small sensitivity to a transverse field when the barrier window is 
located well inside the superconducting electrodes. The possibility to design multijunction 
chips whose each junction has its own magnetic diffraction pattern makes the physics and 
the application of transverse field very attractive and promising. Unfortunately, so far, 
very few experimental works have dealt with transverse field because of lack of theoretical 
understanding. We beheve that this paper will stimulate other groups to fill the gap. 
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